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Abstract 

Problems in lattice gauge models with fermions are discussed. A new 
bosonic Hermitean effective action for lattice QCD with dynamical quarks is 
presented. In distinction of the previous versions Q, it does not include 
constraints and is better suited for Monte-Carlo simulations. 



1 Problems in lattice QCD 



Lattice formulation of QCD seems at present to be the only regular method to study 
nonperturbative effects in this model. Wilson's formulation of lattice Yang-Mills 
theory provides a solid basis for computations in pure gluodynamics. However 
real calculations pose certain requirements on the size and spacing of a lattice. To 
I supress finite volume effects the size of a lattice should be bigger than typical hadron 

diameter, which gives the estimate L > Sfermi. The experience shows that on such 
lattices many physical effects can be described with sufficient accuracy for lattice 
spacing < O.lfermi. 

Modern computer facilities allow to concider the lattices of the size 32^ (record 
size is 64^), which seems to be sufficient for reliable calculations in QCD without 
quarks or in the quenched sector neglecting dynamical quark loops. In this way 
reasonably good results were obtained in calculation of hadron spectrum and form- 
factors, studying QCD at finite temperature and density. 

However, contrary to the comon lore that "all difficulties in gauge theories are 
related to Yang-Mills sector and inclusion of matter fields makes no problem", in 
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lattice QCD the situation is different. Whereas the pure gauge sector is well un- 
derstood, treating of dynamical fermions is far from satisfactory and both principal 
and technical problems are still present. 

Well known problem of fermion spectrum degeneracy on the lattice makes dif- 
ficult to consider chiral fermions as well as to study chiral effects in QCD (For 
detailed discussion and references see e.g. reviews 0, @). Recently a new formu- 
lation of chiral fermion models on the lattice was proposed 0, |]^, which allows 
to remove unwanted fermon doublers preserving chiral invariance of the theory at 
least in the framework of weak coupling expansion. Another approach, pretending 
on nonpertur- 

bative solution of the problem was proposed in papers 0. 

However to estimate a real value of all these proposals extensive nonperturbative 
calculations are needed. Here comes a major technical problem, which is common 
for all models with dynamical fermions, in particular QCD. Calculation of physical 
quantities on the lattice can be done by computing path integrals which are approx- 
imated by corresponding finite sums. But in fermion models one encounters path 
integrals over anticommuting variables which cannot be treated in this manner. 

The action we are interested in is usually quadratic in Fermi fields and the corre- 
sponding integral is Gaussian. It is equal to the determinant of the quadratic form 
in the action, whereas integration over Bose fields produces the inverse determi- 
nant. Hence the most direct way to compute a fermion detreminant is to calculate 
firstly the corresponding bosonic determinant and then to invert it. Alternatively 
one can perform integration over Fermi fields analytically to get a nonlocal bosonic 
effective action, thus avoiding integration over anticommuting variables. This pro- 
cedure includes inversion of huge matrices and is very costly and time consuming. 
However considerable progress in this direction has been done, mainly using hybrid 
Monte-Carlo algorithm H. (For recent review of this and other methods see [M). 



An alternative algorithm was formulated by M.Liischer [0, who proposed 
to calculate a fermion determinant replacing it by an infinite series of boson deter- 
minants. This method is based on the equation of the type 

det(g)2 = jim [det(P„(g2))]-i (1) 

where Pn{s) is a sequence of polynomials such that 

lim P„(s) = < s < 1 (2) 
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By cutting the sequence at some finite n = N one gets an approximate value of 
fermion determinant expressed in terms of integrals over bosonic fields. This idea 
was developped further in papers [T3|, \T^\. It appears however that computational 



efforts needed in this approach are comparable with hybrid Monte Carlo algorithm 
mentioned above. 

Another possibility was discussed recently in our papers , . In our approach 
a D-dimensional lattice fermion determinant was presented as a path integral of 
exponent of D + 1-dimensional constrained bosonic effective action. This represen- 
tation is exact in the limit when the size of extra dimension becomes infinite and 
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lattice spacing zero. For a finite lattice correction terms are present, whose value 
depends on the particular choice of effective bosonic action. In this talk I shall 
use the freedom in the choice of effective action to formulate a new version of the 
bosonization algorithm for lattice QCD which seems to be better suited for Monte 
Carlo simulations, in particular it allows to get rid off the constraint equation. 



2 Bosonic effective action for lattice QCD 

We consider Euclidean QCD on the four-dimensional qubic lattice with spacing 
a. The quark fields are four- dimensional spinors ip^ix), i being a flavour index, 
U^{x) is a gauge field. Let us introduce bosonic fields 4>{x,t) defined on a five- 
dimensional lattice, which have the same spinorial and internal structure as tp. The 
fifth component t is defined on the one- dimensional chain of the length L with the 
lattice spacing b 

L = Nb, 0<n< N (3) 
The effective bosonic action may be written as follows 

N-l 

S = Sw{U) + a' J2{b E {[^"'(C+i(^)0n(x) + - 20:00 - (4) 

X n=0 

-^[0:+i(x)75^0„(x) - h.c.][b]-' + 
+0:(x)D2^„(x) - VLmx)im - tD^,)xix)e-"''^ + h.c.}+ 

Here Sw{U) is the standard Wilson action for lattice Yang-Mills field 0, xi^) are 
four-dimensional lattice bose fields having the same spinorial and internal structure 
as ip. Lattice covariant derivative is denoted by 

D^ij{x) = -[U^{x)tP{x + a^) - ij{x)] (5) 

D = 1/2^,{D; - D,) (6) 

(For simplicity we consider naive fermions, but all the construction is extended in a 
straightforward way to Wilson fermions). 

The effective action (Q) allows to calculate in a standard way any gauge invari- 
ant correlation function. In particular we shall show that the square of the quark 
determinant can be presented as the following path integral 

/ exp{a^ E E '^«(^)(-^+"^)'0«(^)}'^'^'^'0 ~ i ^^v{~S}d(j)*^d(j)ndx*dx (7) 

where free boundary conditions in t are assumed 

(h„ = 0, n<0, n>N (8) 
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For finite b and L this equality has to be corrected by the terms of order 0{h/aY, 
0(exp{— mL}). We consider the case of two fiavours to provide positivity of the 
effective bosonic action. An analogous representation can be written for the modulus 
of quark determinant. 

To prove the equation ( |^ we consider the following integral 



. N-l 

/ exp{a^ ^{6 r'(C+i(a^) exp{-^75Z)6}0„(x) + h.c. - 24>14>^) - (9) 

n=0 



— (x)(m — iD'^^xix) + /i.e.) exp{— m6n}] 
2m 

For small h the exponent in the r.h.s. of eq.( ^ coincides with the effective quark 
action ( Indeed taking into account that the operator 75Z) is bounded, ||75-D|| < 
8a~^ we can approximate exp{— i75D6} by it's Taylor series. Keeping the first three 
terms we get exactly the action ( |). The remaining terms produce corrections of 
order 0(6^/a^), which can be done small by choosing appropriately lattice spacing 
in the extra dimension. 

The operator 75Z} is Hermitean and it's eigenvalues are real. R.h.s. of eq. 
can be written in the basis formed by the eigenvectors of the operator 75-D: 

. N-l 

I = / exp{Y^{b Y: [(r2)(C:iexp{-zD"6}C + h.c. - 20rC) - (10) 

a n=0 

-L^/\(j)^*{m - iD'')x'' + h.c.) exp{-mbn}] 

+ ^X*aXa}}d^Tdrndx"*dx'' 
To calculate the integral (|10D we make the following change of variables: 

< exp{-tD-nb}<P^, exp{zZ}"n6}C* (H) 

Then the integral (p!oD acquires the form 

„ N-l 

I = / expiYib E [b-'iKliK + h.c. - 2</.rC) - (12) 

a n=0 

_^i/2(^^a* exp{-(m - iD'^)bn}{m - 2L'°)x° + h.c] 
2m 

Now the quadratic form in the exponent does not depend on and to calculate 
the integral it is sufficient to find a stationary point of the exponent. For small b 
the sum over n can be replaced by the integral and the stationary equations by the 
following differential ones: 

0*« _ L^/^x*'^im + iD"") exp{-(iD" + m)} = (13) 
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0" - Li/2^"(m - zD") exp{{iD" - m)} = 
0°(L) = 0"(O) = 0, =0"*(O) =0 

The solution of these eq.s is 



(m + iD") 



+ (14) 



l^-i/2 X*"(t-i^) _ ^_i/2 ^exp{-(m + zD")L} 
m + iD" (m + iD°) 



^st 



^1/2 %" exp{-(m - iD'^)t] 



(m - iD°') 



+ (15) 



1^-1/2 X"(^--^) _ ^-i/2 exp{-(m-zD")L}x" 
(m - i/}") (m - iD'^) 

Substituting this solution into eq(|l2D and integrating over t we get 

/ = / exp{- X: ^."^pl^J ^X-rfx" + 0(e-^) (16) 

Integrating over x ^i^nd omitting exponentially small corrections one gets 

I = + = det(-(D)2 + m^) = det(D + mf (17) 

a 

As was discussed above the integral in the r.h.s. of eq.(|^ is a linearized version 
of the integral /, the differences being Oil? Ia}\ At the same time it is easy to show 
that replacing the sum in the eq.( 0) by the integral over t also produces corrections 
of order 0(6^a~^). Therefore we proved that the square of quark determinant indeed 
can be presented as the path integral of the exponent of bosonic effective action (p. 
For a finite lattice the corrections are of the order 

0(6Va') + 0(6-"^^) (18) 

In distinction of the previous version the effective action described here does 
not include constraints, whch makes easier it's Monte-Carlo simulations. At present 
such simulations for simplified lower dimensional model are in progress. 
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